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ABSTRACT 


A method for solving the three-dimensional compressible 
laminar boundary-layer equations for the case of a circular 
cone and a sphere-cone body at an angle of attack is pre- 
sented. The governing equations are modified by a similarity 
type transformation and then transformed into a Crocco-type 
form. The resulting set of equations is solved simultan- 
eously by an iterative method using an implicit finite dif- 
ference scheme by means of an efficient algorithm for equa- 
tions of tridiagonal form. 

The effects of streamline swallowing on a sharp cone 
are included by introducing the true inviscid edge-conditions 
at the distance from the wall equal to the boundary-layer 
thickness. The validity of the approach was established by 
comparison of. the computational results with similar results 
by other methods and with experimental data. 

It was concluded that at sufficiently high Mach number 
and moderate to large angles of attack, the streamline swal- 
lowing effects on a sharp cone result in higher values of 
skin friction and heat transfer as compared with the class- 
ical results for constant entropy. 
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INTRODUCTION 


The Navier Stokes equations of motion which describe 
compressible viscous fluid flow past a body are highly non- 
linear and present formidable obstacles to classical analyt- 
ical or numerical methods of solution except in special cases 
However, since at high Reynolds number the effects of fluid 
viscosity are confined to a thin layer near the body surface, 
many terms in the full equations can be neglected by making 
the boundary-layer approximation. The resulting equations 
can be solved numerically for specified boundary conditions 
at the wall and at the edge of the boundary-layer for pres- 
cribed external inviscid fluid properties. This is possible 
since the boundary- layer equations are parabolic whereas, the 
original Navier Stokes equations are elliptic. Parabolic 
equations are in general easy to solve numerically whereas, 
elliptic equations are not. 

In order to solve the boundary-layer equations, external 
inviscid properties are needed. The usual assumption for the 
external flow at the body surface for supersonic flow is 
that it is of constant entropy corresponding to either the 
oblique shock entropy or stagnation streamline entropy. How- 
ever, as the boundary-layer thickens in the streamwise and 
cross— flow directions in three-dimensional flow, streamlines 
of variable entropy are entrained into the boundary-layer 
and the isentropic condition, especially for a body at an 
angle of attack, does not hold. 

One objective of this study is to develop a method to 
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account for the effect of the variable entropy layer 
near the surface of a body at angle of attack on the solutions 
for laminar three-dimensional viscous flows as compared with 
the classical solutions for viscous flow which have been based 
on constant entropy surface conditions. It is interesting to 
note that despite a large body of literature related to three- 
dimensional boundary layers, this particular area has remained 
relatively unexplored. 

The literature reviews by Cooke and Hall (ref. 1), Mager 
(ref. 2) and Stewartson (ref. 3) reflect the major extent of 
work done in the area of three-dimensional boundary-layers. 

A case which has received considerable attention, because of 
simplifications due to spanwise derivatives being zero, is 
conical flow. The majority of the work done on viscous con- 
ical flows has been concerned with approximate methods of 
solution which are usually valid for small cross-flow and 
small angle of attack. Only recently the restriction of small 
cross-flow has been removed in work for sharp circular cones 
by Boericke (ref. 4 ),. McGowan and Davis (ref. 5), and others. 

In conical flows the boundary-layer solution is actually two- 
dimensional since similarity exists in the streamwise direction. 
The flow near the windward side creates no problems and is 
dominated by inviscid conical flow with a favorable pressure 
gradient. The region near the leeward plane at high angles 
of attack may be strongly influenced by inviscid viscous inter- 
actions. There is flow turning toward the windward plane of 
symmetry which may give rise to imbedded shocks if the local 



Mach number in the cross-flow plane is supersonic. This may 
cause separation. The presence of an adverse pressure grad- 
ient can cause the cross-flow velocity to reverse and cause 
the flow to recirculate. Typically, the flow in this region 
even at moderate angles of attack, exhibits a bump with stream- 
lines diverging as they approach the leeside of the body, 

(Fig. 1.1). 

A cone at angle of attack is characterized by the fact 
that the streamlines around the periphery of the cone are 
passing through a shock of varying strength and posses diff- 
erent entropies. The surface streamlines pass through the 
windward meridian plane and have the highest value of entropy. 
The streamlines away from the windward plane pass through a weak 
er shock and posses lower values of entropy and enter the 
flow field nearly meridionally. There is a layer near the 
surface in which the entropy varies rapidly with the distance 
from the surface, from the maximum value, to a lower value con- x 
sistent with a weaker shock. Thus, although the entropy at 
the surface is constant, the entropy of the streamlines en- 
trained by the boundary-layer is of varying magnitude. The 
question of the magnitude of the effect of the variable entropy 
on the viscous boundary-layer flow has been of considerable 
concern and importance. This problem also occurs on a blunted 
body at zero angle of attack with a bow shock which is normal 
for the stagnation streamline wetting the surface. The 
streamlines passing through an oblique shock have lower entropy 
as they are "swallowed" by the boundary-layer downstream. 
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(Fig. 1.2). The same phenomenon occurs on a sharp cone at 
angle of attack (Fig. 1.3) where the streamlines away from 
windward plane enter the boundary-layer with lower entropy 
and create a thin vorticity layer. Of course, outer edge 
boundary conditions accounting for streamline swallowing will 
yield results different from the classical approach using edge 
conditions of constant state as was concluded by Ferri (ref. 

6) and Mayne and Adams (ref. 7). A more proper approach to 
account for swallowing of a streamline as it enters the boun- 
dary-layer, would be to use the edge conditions consistent 
with the flow conditions at the distance from the surface 
equal to the boundary-layer thickness. 

In order to establish the extent of the influence of 
the streamline swallowing upon the boundary-layer, a method 
for accounting for the effects of vorticity swallowing has 
been developed. The objective of the technique presented 
here is to determine the effects of the variable entropy edge 
conditions on solutions to the three-dimensional boundary- 
layer equations for supersonic flow over a right circular cone 
at angle of attack. The second part of this work is devoted 
to the solution of viscous laminar flows over blunted cones. 

The governing boundary-layer equations will be first 
modified by a similarity type transformation and then trans- 
formed into a Crocco-type form. The resulting equations will 
then be solved simultaneously by an iterative method using an 
implicit finite difference scheme similar to one developed by 
Dwyer (ref. 8) with modifications by Kraus (ref. 9). 
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2. THREE-DIMENSIONAL COMPRESSIBLE LAMINAR BOUNDARY- LAYER 
EQUATIONS* 

An implicit finite-difference technique originally de- 
veloped by McGowan and Davis (ref. 5) will be employed for 
solving the three-dimensional compressible laminar boundary- 
layer equations over sharp and blunt right circular cones. 
The governing equations for a generalized body-oriented co- 
ordinate system to be used in the analysis (see Fig. 2.1) 
are as follows: 


£ - momentum 


u_ au + v_ 3 u + w_ 3u + uv 3h l _ _2L_ 

h ■ 3£ h 3n h 3? h h 3n " h h 35 

h l n 2 n 3 n l n 2 n l n 2 



^ ?h 3 




( 2 . 01 ) 


n - momentum 


u_ a v , v_ av w_ a_v 

r 3 C H • 3 n r 3 C 

h l h 2 h 3 



h l h 2 





Sh 2 

3 "? 



( 2 . 02 ) 


* Detailed presentation of the equations and transformations 
discussed in this section is given by McGowan and Davis 
(ref. 5) and Der (ref. 10). 
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energy 


u at v_ at. w at 




3 n C 35 


V Lh i 


u ap v ap 


E 3 5 c 3n 
2 


(2.03) 


+ 


V S 3 


fau ] 2 [ av 1 

La c J L3 5J 


- 2 


C ph-. 3C 
P 3 


continuity 


k at 
h 

3 


ff <h 2 V 5) + < h iV 5) + |j < K i fi 2 pw) 


= 0 , (2.04) 


state 


p = Rpt 


(2.05) 


and 


viscosity (Sutherland) 


3/2 




fc r + C 


t + C 


-V 


(2.06) 


where C = 198.6 R for air. 

The ^-coordinate is chosen to be the arc length along 

r 

i. _ _ _ __ 

the body; therefore, h^ will be a constant, L. and are 
chosen to be arbitrary functions of the coordinates. £ , nz 
and c are dimensionless, consequently the scale factors h^. 
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h 2 and have the dimensions of length. Similarity- type 

transformations for the normal coordinate and velocity are 
defined as follows: 



(2.07) 


and 



(2.08) 


where 5 * is an arbitrary function of 5 and n and will be 
chosen such that similarity conditions are satisfied. For 
convenience the following definitions are made: 


F = — , G = — , e=£-,H= 2L and £ = | . (2.09-2.13) 

e v e e e 

The equation of state in the dimensionless variables is 
written as: 


0 = V‘ • 


(2.14) 



- 8 - 


The scale factors h 2 and are defined such that the 
equations will be in a Crocco-type form. The dimensionless 
scale factor h 2 is a function of £ and n? that is 

h 2 = h 2 U/n) • " (2.15) 

The normal coordinate 5 is defined as 

C = /I - F (2.16) 


and finally we choose the shear parameter <j> to replace F 
as a dependent variables as 


. £ 1 3F _ £ 

* h 3 k “ " KJ 


(2.17) 


Equations (2.09) to (2.17) are next substituted into the 
governing equations (equations (2.01) to (2.04)). These 
variables are similar to those used by Der (ref. 10) . The 
main advantage of these variables is that in the numerical 
calculations the integration is always carried out between 
C = 0 (f reestream) to 5 = 1 (wall boundary) . Thus use of a 
specified finite region in the ^-direction can save large 
amounts of computer time over other variables which extend 
over an infinite region. This is important in three dimen- 
sional boundary-layer calculations where computing times 
can become excessive. 

The governing equations (equations (2.01) to (2.04)) can 
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then bo expressed in terms of Crocco-type variables as 
follows: 


A. (1-s 2 ) 2 + A c (l-£ 2 ) G + A, ^ - A, G 2 = 


" A 7 e " f It (<k) 


(2.18) 


(1-C 2 ) 


3G 

35 


a 2 g 


A 1O 0 + t fc 


3G 

3n 


U 


- A. 


3G, 

3? J 


4>H 

£ 


a^G 

3? 


A 8 (1 " C< 


) G + 


A g G 


(2.19) 


At (l-S 2 ) 14 + A 2 G |4 " A, ^ |4 - [A^-A^] (1-c 2 ) 


35 


3 • £ 3 £ 


'll 18‘ 


(2.20) 

+ !A 12 -A 19 ,G0 + [A 13 c 2 + A 14 (|f) 2 ] 4 + h T k <* H’ 


and 


V 1 -? 2 * It <t> + A 2 k (l f’ - a 3 If + (4A 4 + A 1S> (1 -c 2 > 


( 2 . 21 ) 


+ (2A 17 + A 16 ) ^ 0. 


The coefficients A^, A 2 , ... , A 19 are functions of the 
inviscid quantities and body geometry and are defined in 
Appendix A by equations Al to A19. In order to reduce the 
system to three governing equations, equation (2.18) is solved 
for H as follows: 


l-e- 
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AgG 2 !. A 5 G(1-C 2 )2- A 4 (l-C 2 ) 2 «. A^QSL 

A 3 C$ A 3 C* ; A 3 5«j> Ajc? 


l a_ 

A 3 C 3? 


(<K) 


( 2 . 22 ) 


Substituting equation (2.22) into equations (2.19) to (2.21) 
leads to the final form of the three governing equations for 
the dependent variables G, 0, and s 




lv 3G 

? 3? 


+ B 2 G + B 3 + b 4 


3G 

3C 


+ B r 


3G 

3n 



(2.23) 


lie 

3C 2 


+ t (1 Pr) li-l£.+ B, ] + B^0 + B fi + 

d> 3 C C 6 3c 7 8 


B 


30 
9 35 


+ B 10 = 0 - 


(2.24) 


and 


3 2 <)> . ( 1 . B 11 i 3* . <J> . B 12 . B 13 3 <J> B 14 3 <J> 

77 « 7~ 7 * 77?? + *?"< 


= 0 


(2.25) 


where B^, B 2 , ... , B^ 4 are defined in Appendix A by equations 
A20 to A33. 

The viscosity function appearing in the B^ coefficients 
(see equations (A20) to (A33)) is given by Sutherland's law as 


l 


1 +• c /t r 
9 + C/t e 



1/2 


£ 

t 


(2.26) 
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The quantity a?,/a t (see equation (A31) ) is obtained by dif- 
ferentiating equation (2.26) with respect to t; that is, 


') _ 3 9, 'd 6 • 

5^ D0 jt 

The boundary conditions on G and e at the body surface 

are : 

G = 0 (2.27) 


t 

6 = — , at 5 = 1 (2.28) 

e 

or 


M 

3? 


f U *r\) 


(2.29) 


where f(£, n ) is related to some specified body surface heat 
transfer, and t^ is some specified body surface temperature. 

At the outer edge of the boundary- layer the boundary 
conditions on G and q are: 


G = 1 




(2.30) 


and 


f at c = 0 


0 = 1 .. 


(2.31) 


J 

The boundary conditions on <j> are more complicated. 
Applying L' Hospital's rule to (2.17) evaluated at the outer 
edge of the boundary-layer yields: 




3 2 F 
3? 2 


(2.32) 


C = 0 
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a F ... 

The quantity - — 5 - must vanish at the mviscid edge; therefore, 
* H 

at £ = 0 

4> I =0 (2.33) 

lo 

The boundary condition on 41 at the wall can be derived from 
(2.18) applied at the wall. Setting 5 = 1 and using the no 
slip conditions result in: 

A 3 H * A 7 *7 - ♦>„ (2 ‘ 34) 


The nondimens ional heat- transfer rate can be shown to be: 


Ul = ( - 

a «J.w- 


P « 


( -v ) 

Y-l 


t 

e U 1®. ) 

P r g /Re w 


(2.35) 


Also the skin friction coefficients can be shown to be: 



y 3U 





4 Vw u e 

{*/Re 


and 


(2.36) 



2 Vw v e 3 G\ 

S*/Re 35; ' 

^ w 


(2.37) 



- 13 - 


Finally, the physical normal distance can be expressed as: 


1 

— , 9 * 

s 3 “ “ / J « ed ? 

c 



( 2 . 38 ) 


The governing equations and boundary conditions are 
now in a convenient form for solution using numerical methods . 
The next section will deal with developing an implicit finite 
difference method of solution. 
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3. NUMERICAL METHOD OF SOLUTION 

The governing equations presented in Section 2 will no w 
be reduced to finite-difference form in order that they may 
be solved numerically. The complete solution to be obtained 
on a digital computer will be comprised of three separate 
cases: a general scheme, a scheme in which similarity exists 

in one of the surface variables, and a scheme in which simi- 
larity exists in both surface variables. Finally, the com- 
plete computation scheme will be explained as a combination 
of the three separate cases. This development is identical 
to that presented in reference 5. 

The governing equations are solved using a step-by-step 
finite-difference method similar to that suggested by Dwyer 

(ref. 8) with modifications by Krause (ref. 9). The method 

2 

has accuracy of order A , where a is either A £ or An> and is 
stable for negative transverse velocity if proper step sizes, 
are chosen. A schematic of the difference model is presented 
in Figure 3. a. It should be noted that the equations are 
written at the point (0,n) and solved for the values of the 
quantities <j>, e and G at the point (2, n) . Thus, the quantities 

2 

(1 -f } ' (t 7 } ' (If) ' ' and “n n' where the 

B C 2 0 , n 85 0,n 0 ,n 8t1 0,n °' n 

subscripts (0,n) refer to the point (0,n), can be expressed 
as : 
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o) _ w 2,n + “3,n 0 , . 2 

0 , n ~2 + 0{A >' 


f — ^ \ +( 5jE±L \ 

( iuA _V a ^2,n V 3 ^3,n . 

^i 0 ,n' 2 


0 (A 2 ) , 


and 




/ 2 , 

(M) 

' 3 t / . 


+(M) 


^ y 2,n Vg g 


^ + 0 (a 2 ) , 


where 


ltt) i>n v i '% t i '°' 1 - °^ 2 ' 


and 


(4). 

v 3c /i,n 


C with 


(l) • > 2 (jJ • I (0 • <y »*i 

£±i hj* i,n z l + 0 (Ai; 2 ), 

AC J 


( 


3m \ 

9 ^0,n 


^2 , n u 3,n + „ ,.2, 


AC 


0(A) , 


and 


( 

V 9 Vo,n 


(<0 o • — OJ t _ ) + (tO. . “ (II, n ) O 

- 2 ' n L'» Ln_ + 0 (a 2 } 

2An 


(3.01) 

(3.02) 

(3.03) 

(3.04) 
i = 2,3 

(3.05) 

(3.06) 


(3.07) 
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The governing equations (2.23) through (2.25) can be 
written in the following general form: 


(MO + a l (ff) + a 2 u n + a 3 + a 4 (ft) + a 5 (ft) 

' 3 £ 'n n' n n n v s/ n n> M/ n 


= 0 


(3.08) 

Substituting equations (3.01) through (3.07) and re- 
arranging yields the following expressions for (3.08): 


A (l) A» 1 I 5 A I C 1 1 \ A . *1 — * . 

n 2,n-l n 2,n n“2,n+l n 


(3.09) 


where 


"1 1 
= - (— 


a l 


n 


n 


2 ' 2 
A? 


2 A £ 


) , 


(3.10) 


a, 


= (- 


n 


n 


» 

a 5 


n 


A ? 


A6 2An 


-) /■ 
n ' 


(3.11) 


a l 

C = - (— + — 

n 2 1 2 2 a ? 


) , 


(3.12) 


and 


D = 
n 


Clo .+ (~ — + 


n 


2 r 4 a?' W 3,n+1 


n 2 a? 


(3.13) 
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ct 


+ ( — — 
2 a c 


n, 


Th7 ) 0) 3,n-l 


I 

u 5 

n 

2a n “l,n 


a 2 


+ (- 


n 


n 




AC 


AC 


2ah / w 3,n n 2An w 4 ,n 


n 


a) , 


In the special case for which there is similarity in the 

C-direction for example at the tip of a sharp cone or the 

nose of a blunted cone, a Crank-Nicolson finite-difference 

scheme can be used in the cross-flow direction with accuracy 
2 

of order An .In this scheme, the equations are written at 
the 0-point and solved for the value of the quantities e, <J>, 
and G at the 2-point. Thus, yields the following relations 
(Fig. 3.1b) for the value and partial derivatives of a vari- 


able : 


_ w 2,n + w l,n , 2, 

w 0,n j ' + 0(An )r 


(3.14) 



(3.15) 


(3.16) 
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where 


j = — L ' n+ 1 — — + 0(Ac 2 ) 
! '4,n 24 « 


and 


(3.17) 


V with i = 1,2 


(l co \ _ W i, 

' 3^ 2 A,n 


,. -201. +o). . „ 

s+i il“ + 0{fi ? 2 ) 


Ac' 


and 


) 


(3.18) 


' 3 n ' 0 , n 


= - 2 . ^ kdl + o (An 2 ) . 

An 


- 0) . 


(3.19) 


Similarly, as before, the governing equations can be 

written in the form of equation (3.08) with = 0 ; upon ■ 

n 

substituting the finite-difference form of the partial de- 
rivatives as specified in equations (3.14) through (3.19) 
into (3.08) a relation of the form of equation (3.09) is 
obtained. In this case. 


A = (- 
n 


n 


2 Ac 


4 A c 


) , 


(3.20) 


B = ( 9 

n 2 


n 


1 + 


AC 


2 An 


2 


(3.21) 
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C = { 
n 


n 


2a c ' 


4ac 


) , 


(3.22) 


and 


a ■ 


D = - 
n 


a 3 


' + (-4 + 


n, 


n 


2ac 


4ac’ } “l,n+l 


(3.23) 


a- 


n, 


i 

a 2 


L 2a c‘ 


4ac ' w l,n-l 


+ (- 


n 


AC 


n, 


2 An ' “l,n 


Equation (3.08) can be evaluated in a similar way for the case 

I 

similarity in the n-direction with a 5 = 0. 

n 

Finally, an iterative finite-difference scheme is used 

for the case of similarity in both the ^-coordinate and the 

n-coordinate. In this case the governing equations reduce 

to ordinary differential equations with the equations being 

written at the 2-point and the values generated at the 2- 

. 2 

point; subsequently, the accuracy is of order Ac . The par- 
tial derivatives, which are now actually total derivatives, 
become : 



2,n+l 2 * 11-1 

2A C 


+ 0 (AC; 


and 


(3.24) 


( 


T 2 - 
d W ^ 

9C 2 ' 2 , n 


u>~ . — 2u)_ + 

2 , n+1 2 , n 


U) 


2,n - ~ + 0 (Ac 2 ) 


Ac‘ 


(3.25) 
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Again, the governing equations can be written in the form 

1 I 

of equation (3.08) with = a 5 " Upon substituting the 

n n 

finite-difference form of the partial derivatives as expressed 
by equations (3.24) and (3.25) into equation (3.08) equa- 
tion (3.09) results, where: 

I 

tt l 

A n = - 3sf> ' (3 ‘ 26 > 

A 5 


B n = ( “^2 + a 2 )r ' (3.27) 

A? n 


i a l 

C n “ -h + ' ' (3 - 28 > 

A £ ^ : 


and D n = - o 3 • (3.29) 

Thus, one can apply equation (3.09) at N-2 evenly spaced 
points in the normal direction from point 2 , one step away 
from the outer edge, to point N-l one step away from the wall. 

The general solution to equation (3.09), as given by 
Richtmeyer (ref. 11) for example, is of the form: 


03 


2 , n 


E n w 2 ,n+l 


1 

+ F 

n 


(3.30) 


If one applies the outer-edge boundary conditions at 
5 = 0 and requires that equation (3.30) holds independently 
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of step size, the resulting conditions are: 


and F 1 = 0,1,1 respectively for <J> , 6, and G. 

Then the following recursion relations can be derived: 


(3.31) 


(3.32) 


With the known values of the variables at the wall and 

I I 

known E and F , one can calculate the values of . to 

n n ' 2 ,N— 1 

w 2 , 2 * 

*> 

t 

The boundary conditions at the wall for 0 and G are ^ 

e 

and 0, respectively. However, the boundary condition at the 
wall on $ from (2.34) has the following form: 


and 




— C 


E = 
n 


n 


A E , + B 
n n-1 n 


• ii 

, D - A F , 

' _ n n n-1 

F r ~1 -r 

A E , + B 
n n-1 n 



{ A 3 H 2,N 


A 7 £ 2?N 2 ,N _ 
^2 ,N 


^2 ,N* * 


(3.33) 


With a 3-point difference form for the above partial 
derivatives at the wall, one obtains: 


[?± 


/ 2 , N 


2 (3<j> 2,N 


4 ^2 , N-1 + 


<f) 2,N-2 ) 


0(A^ 2 ) 


(3.34) 
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Now using equation (3.09) for $ at the point N-l gives 
and additional relation: 

^SsJ-1^2 ,N-2 + B N-1 < **2 ,N-1 + C N-1^2 ,N = °N-1 (3.35) 

The final expression for this boundary condition is determined 
by applying equation (3.30) for <p at N-l: 


^2,N-1 E N-1^2 ,N + F N-1 


(3.36) 


Applying equations (3.34) through (3.36) and evaluating 
for <J>2 n ' one arr; *- ves a *- a Quadratic equation having solution 


-b’ - 


f I 

4a c 


2 ,N 


2a 


(3.37) 


with 


• C N-1 ' ' B N-1 

a = _ + 4E n . 1 + -r-i 


- 2AC - 3 . , 


V: 


^-1 


(3.38) 


B. 


b = 


4f n-i + f n-i 


N-l D N-1 


^-i ^Sj-: 


A 3 2ACH 2,N 


and 


(3.39) 


C ~ A 7 £ 2,N 6 2,N 2A? * 


(3.40) 
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It should be remembered that the coefficients A in 

n 

equation (3.09) are functions of the dependent variables 
0, and G and are unknown a priori. An initial guess for these 
unknown variables is required and the result is iterated 
upon until the values of these unknown variables are relaxed. 

It has been found that the 1/ 4> term in equation (2.25) 
and subsequently in equation (3.08) caused the iteration 
scheme to diverge. Therefore, a "quasi-linearization" of 
this particular term was introduced as follows: 

1/4. = 2/ 4> g - 4>/4> g 2 (3.41) 

where 4> G is the value of the shear parameter at the previous 
iteration; with this modification the solution coverged 
rapidly. • . 

3.1 General Solution Technique 

The system of non-linear second order parabolic equations 
developed in Section 2 are solved simultaneously by an 
iterative process for a sharp circular cone and for a spheri- 
cally blunted cone. The governing equations ( 2 . 23 ) to (2.25) 
are reduced by replacing the partial derivatives by finite 
difference expressions, to algebraic equations having the 
tridiagonal form of equation (3.08). The coefficients in 
equation (3.08) are functions of the dependent variables 0, 

and G which are to be determined and must be represented 
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by their assumed initial values and are subsequently updated 
during the iteration process. At the tip of a sharp cone or 
at the stagnation point of a sphere cone (which correspond 
to ( = 0 and n = 0) similarity in the £ and n directions 
exist. The mathematical solution at this point is considered 
as a limiting case as ( + 0 and n -*• 0. The governing equations 
at this point reduce to ordinary differential equations in 
£ only. The solution to this two-point boundary-value prob- 
lem is a starting solution. For grid points on the surface 
corresponding to £ = 0 and n > 0, similarity in £ exists (there 
is no £-dpendance in the equations) and the solutions are 
obtained using the starting solution and marching circumfer- 
entially around the body in the n-direction. Thus, the sol- 
ution in the (n,£) plane is reduced to a two-dimensional 
problem. In the plane (£,£), with n = 0 and £ > 0, similarity 
in r\ exists and two-dimensional solutions in the windward ' 
meridian plane are obtained. For points outside of the £ = 0 
and n = 0 planes the dependance on £ , n and £ exists and the 
solutions are obtained by using the solution at the windward 
plane and marching in planes around the body in the n-direction 
and stepping forward in £-direction (Fig. 3.2). At each 
(£,n) grid point on the surface, iteration is applied on the 
values of 8 , <j> and G in the normal direction until the slope 
of the shear parameter in the normal direction is relaxed. 
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4 . GOLUT [OWS FOR SHARI* CIRCULAR CONKS AT ANGLE OF ATTACK 
WITH STREAMLINE SWALLOWING 

The recent contribution in the field of three-dimension- 
al boundary-layer studies by McGowan and Davis (ref. 5) has 
extended the numerical analysis of three-dimensional laminar 
flow to sharp circular or elliptical cones at angle of 
attack. The work by Boericke (ref. 4) treats laminar flow 
over circular cones at angle of attack using Moore's trans- 
formation and numerical as well as experimental surface 
pressure distributions. Quite recently the work by Mayne 
(ref. 12 ) presents an analysis of laminar boundary- layers 
on right circular cones at angle of attack, including the 
effects of streamline entropy swallowing with the assumption 
of a negligible effect of streamline swallowing in the 
windward plane (a similarity solution at £ = n = 0 has been 
used to apply to the region 5 > 0 , n = 0 ). 

In the present work a method for solving for flows over 
sharp cones at angle of attack with streamline entropy swal- 
lowing is developed. The details of this method will be 
presented in the following sections. 

4.1 The Geometry of the Problem 

In spherical coordinates, an infinitesimal arc length 
is given by (Fig. 4.1): 

d£^ = dr^ + r^ (d$^ + sin^de^). 


(4.01) 
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l.n the £, n , c body-oriented coordinate system, as shown 
in Figure 4.2, an infinitesimal arc length along the surface 

(4.02) 


(4.03 a-f) 


4 . 2 Conditions at the Outer Edge of the Boundary Layer 

The usual procedure of matching the viscous and inviscid 
flow fields at the edge of the boundary layer assumes the in- 
viscid value of the pressure at the body surface and the en- 
tropy value consistent with the oblique shock or stagnation 
streamline value. It has been indicated, however, that there 
is a need to develop solution techniques which include shock 
layer vorticity (ref. 13). In the current scheme the initial 
step in determining the inviscid edge properties is to obtain 
a boundary layer solution for inviscid values at the surface. 
Having obtained the magnitude of the boundary-layer thickness, 
new inviscid values are obtained by interpolating the inviscid 
properties for cones (ref. 14) around the body at a distance 


of a sharp cone is 


2 2 2 2 2 
djT = h r z dE + h 2 d n 


where one can identify 


r = 5 and ~ K 


h^ = 1, dn = sinc() dG 


3h~ 3h_ 

= 1 and =■ = 0 


3C 


an 
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from the wall equal to the boundary layer thickness. With 
these inviscid properties, a new boundary-layer solution is 
then obtained at the next downstream station. 

The scheme adopted for determining the inviscid properties 
is based on the invariance of total enthalphy in the inviscid 
part of the shock layer which is true if the total enthalpy 
in the freestream is constant. Due to the presence of vor- 
ticity in the inviscid flow, the velocity and temperature 
gradients are no longer zero at the boundary layer edge. How- 
ever, it is assumed that the net energy flux between adjacent 
fluid laminae in the edge region, particularly on sharp bodies 
becomes negligible. 


4 . 3 Non-Isen tropic Conditions 

The entropy is defined by the expression 


ds - <£i) 

T rev 


Using the thermodynamic relation 


ds = 

m 


one can show for a perfect gas with constant specific heats 


with 
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pv = RT , 



that the change of entropy is related to pressure and density 
by the expression 


S - S c - — C £n 
ref v 


P/P 


ref 


(p/0 ref ) Y 


Therefore the parameter 



(4.04) 


is a measure of the change of entropy and can be used to re- 
late the value of entropy to pressure and density. 


4 . 4 Boundary Layer Thickness 

The physical coordinate normal to the surface of the 
body is defined by equation (2.38) as 



1 

“ " / J e ^ 

c 


L 

/Re 


which results in 


S /Re 

: = - / f edc = Y(c) 

L F, 9 


(4.05) 


(4.06) 
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If we define the thickness of the boundary layer at the point 

where the local tangential velocity u in the streamwise £- 

direction reaches 90% of the edge value or where F = — = .9, 

e 

we obtain the boundary layer thickness 6 = __ Q : 


6 /Re 

L rc 


= y(N) 


F = 


.9 



where £ was replaced by /g in order to insure similarity at 
the tip (see section 4.9). With Reynolds number defined by 


co ■» 

we obtain an expression for the boundary layer thickness 
- - V? Y „ = q l 

6 = S-, = - — (4.07) 

^ /Re 

4 . 5 Location of the Edge of the Boundary-Layer 

The location of the edge of the boundary-layer in the 
shock layer is specified in terms of the \ parameter which 
identifies rays in the conical flow (Fig. 4.3): 

l = * 

b 

The expression for the \ parameter when evaluated in terms 
of the boundary- layer thickness 5 and body coordinate 


takes the form 



-30- 


2 

6/cos d> 

Z = — r 2 (4.08) 

(CL - 6 tan <j> ) (tan $ [<p) tan <p ) 

C 5 C 

The shock angle <j> is interpolated by the interpolation rou- 

O 

tine as <p = <f> (n) and the 6 value is known from the previous 
s s 

C station. 


4.6 Inviscid Data 

The inviscid data used as an input to solve the boundary- 
layer equations is obtained from Jones (ref. 14). In order 
to obtain a consistent set of inviscid data for the case of 
entropy swallowing, the input values of v g , w q , p e , and p e 
at discrete points around the periphery of a body, and 
Bernoulli's equation 



2 2 
+ w + u + 
e e 




(4.09) 


are used. 

In the above equation the expression for v fi is found 
from the c -m o ment u in equation 


9u cu au „ 9p 

v = — — + — — - — — + — § — 

e 8ri v e v e p e 


(4.10) 


9 U 


. . 0 

and a 5-pomt difference expression for is used, i.e. 

9 n 


SU G _ 8 (U i + l~ U -i-l ) 




9 n 


12 a n 


(4.11) 
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The follov/ing two expressions require special attention 
in the plane of symmetry where v g = 0 : 


l_!^e _ 

v e 3 5 v e 3 5 


(4.12) 


The first expression is obtained by a limit process and the 
second one is obtained from the 5-momentum equation 


5 3u e _ 1 


v e 


(v 
u e e 


JL_ 1_ !£e 

3n P e v e 35 


(4.13) 


With equation (4.10), equation (4.09) becomes 


(!M 2 + u 2 + u 2 U !V)\ (g_ u !^e f2_, 

\ 3n / e e \ v e 35 / VP e v e 9 «./ e 3 n V v e 9 C / 


3 U ^\2 /f . 3P_\ 2 3u„ 3u 


»n \ P e v e 35 y e \.v e 35 / \P e v e 35 ) 


, £) + w 2 

35 ✓ ■ e 


+ 2 




r — 1 1 

V KV- ho ) = 


(4.14) 


which is of the form 


(v u 


3“ e 3u e 3p e 3p e\ 
e ' V e ' w e' p e' ' 35 ' 3 n ' 35 ' 3 


n ) " ° 
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Introducing equation (4.11) and evaluating the equation (4.14) 
at N points/ a set of N simultaneous nonlinear equations in 


u = (u^ u 2 , ... , u n ) 


is obtained and is solved by the Newton-Rapson iteration tech- 
nique using the Gauss- Jordan method. 

Defining 


F. . (u) 
xj 


3F i (.u) 

9U . 


= <j> (u) 


(4.15) 


with the starting velocity vector 


u 0 ~ (u 10' u 20 ' *** ' u n0 } 


the solution vector for velocity is 


- k+1 - K , t K 

u = u +6 


(4.16) 


__ K 

where 6 is the solution vector for the set of simultaneous 
linear equations 


4> (u) 6. = - F i (u) 


(4.17) 


and is obtained by the Gauss- Jordan method. 
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With 

6 j = - AU j (4.18) 

the solution vector for velocity is determined for the k-th 
iteration 

u. k+1 = u. K - &Uj K . (4.19) 

In order to compute the values of u g which are consistent 
with the entire set of the inviscid data, the following in- 
viscid values have to be given or generated: 

3p_ 3 u 3u 

u e' v e' 'V K “ jf' ^ 5T' 

where u^ can be used as initial values. In evaluating the 

inviscid values, the symmetry relations for the windward and 

leeward planes are utilized. All interpolations are made by 

the cubic spline method (Appendix B) . 

In the present analysis, including the effect of the 

entropy swallowing, the edge conditions are determined at 

the normal distance <5 from the wall and they vary from station 

SO 

to station in the £ direction giving rise to gradients —• , 

3v e 

where Q = u , v , p , and t . The derivative -r is obtained 
e e' *e e 3£ 

from the n-momentum equation (2.02) as 
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ac 




e + 1 


3P, 


V V \ 
e e \ 

3 n / 


Z SP e u e 3n 


( 4 . 20 ) 


3t e 

and the — - — derivative is obtained from the energy equation 

d 


(2.03) as 


3t e 

1 

" r± 

1 (u 8p e + v e^e\ 

v 3t 1 
e e 

3 Z 

u 

e 

Y 

p e V e H Z 3n/ 

Z 3n 


( 4 . 21 ) 

8t e ; 

In order to obtain - xn this expression, the equation 

3 n 



is used as follows: 





+ 



3K 
3 n 


(4.22) ~ 


For the case of constant entropy with zero values of deriva- 

9 0 

tives in the quantities , the inviscid E; -momentum equation 

0 c, 

applied at the surface yields: . 



(4.23) 


and using the fact that the entropy is constant on the sur- 
face of a cone: 
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(4.24) 


One can then generate the entire set of inviscid data by 
the method described previously. 

The value of the inviscid, outer edge temperature is 
obtained from (4.18) along with the perfect gas law (2.05) to 


t e = (KpJ-V^ , 

and therefore 

3t , . t 3p 

e _ (y — 1) e e 

a n y P e 3n 


(4.25) 


(4.26) 


From the inviscid n-momentum equation, one can get an ex- 
pression for the pressure derivative: 



, 9 1 * * * V e 
v e ( 7fT 


+ u e> 


(4.27) 


This form of the pressure derivative is convenient since the 
1 


value of 


v e 3n 


P e s P e 

is necessary in (A10) and although both 


a n 


3 P, 


1 ^ e 

and v g go to zero at certain points, the quantity — 3 - - 
does not, for example at the windward streamline. 

The inviscid flow properties are determined by the 

previously described method wich reduces to a solution of 

a set of non-linear simultaneous equations for u e> In solv- 
ing this system the following two cases are considered: 
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1) for isentropic conditions the surface entropy is 
assumed to be constant for the entire shock layer, K - p/p Y 

constant 

2) for entropy swallowing the entropy varies with n, 

5 and c, therefore, for each C station the value of K at the 

edge of the boundary layer ' is K = K(n) = — . Since both 

P (n) Y 

pressure and density vary in the normal and circumferential 
directions, p = p(n,0 and p = (n,s) are obtained by a double 
interpolation routine using cubic splines as p(n) = p(n><5) 
and p(n) = p(n,6). With the p(n) and p (n) values known, the 
velocity u e is obtained by the method described in this 
section. In order to make the inviscid values generated by 
Jones (ref. 14) consistent with this program, certain con- 
version factors have to be introduced: 

p = — = 0 ^, where subscript J denotes data by Jones 

p j 

oo 

and 

p = E_ = Pj (l m 2 sin 2 <J> c ) 

4 - V Limitations of the Inviscid Flow Methods 

The boundary layer equations can be solved once the 
inviscid edge conditions are prescribed. The computation 
of the inviscid properties is complicated by the existence 
of a strong entropy gradient and the occurrence of a singu- 
larity at the leeward symmetry plane. 
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The flow around a cone at an angle of attack, which is 
of interest in this work, is rather complex. One can dis- 
tinguish the region on the lower surface of the cone which 
can be represented by elliptic equations. In this region 
the effect of the boundary layer is present but one can 
assume that the results of the inviscid theory are valid. 

The flow above the body is represented by hyperbolic equations 
and can be solved by the method of characteristics reaching 
in some cases as far as leeward plane. However, the flow 
symmetry is not satisfied in this region for high angles of 
attack which indicates in general the possible presence of 
internal shocks (ref. 27). The region near the leeward plane 
with separated flow and internal imbedded shocks is dominated 
by strong viscous effects and should be handled by equations 
which include viscous effects even outside of the boundary 
layer (ref. 15). The possibility of the presence of internal 
shocks can be inferred from the fact that the flow Mach number 
in the cross-plane reaches values greater than unity on the 
upper surface (Fig. 4.4). It can be concluded from the scatter 
of the inviscid data that the computed inviscid data are not 
accurate for these conditions and the computed viscous para- 
meters may not be precise. 

Current techniques for obtaining the inviscid flow 
properties include perturbation methods (refs. 16-18) and 
numerical methods (refs. 19-22). The perturbation methods 
rely on first order perturbation solutions for the velocity 
components and lead to considerable inaccuracy in velocities 


I 
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and pressure at large Mach numbers and at moderate to large 
angles of attack. The numerical methods for computing the 
inviscid conical flow are applicable at relative angles of 
attack a/<j> c up to unity. The computed pressure distribution 
is generally fairly accurate, however, the computed values 
of velocity components u and v are not. A more consistent 
set of inviscid data is possible by using the computed or 
experimental pressure distribution and entropy value at the 
windward plane and computing v, p , and u from cross-flow 
momentum, conservation of total enthalphy and constancy of 
entropy at the surface (ref. 4). A comparison of the pressure 
distributions computed by various methods (ref. 4) including 
the values used in this analysis is given in Figure 4.5. 

The pressure by Moretti (ref. 19) and by Jones (ref. 14) 
give the best agreement with the test data by Tracy (ref. 23) . 
Since the present method accounts for the effects of entropy 

v> 

swallowing and the pressure distribution across the entire 
shock layer is required, the readily available inviscid values 
by Jones were used. It should be noted that this pressure 
is a fairly good approximation to experimental pressure over 
most of the body, but may become inadequate in the leeward 
region where the strong interaction with the viscous flow is 
not accounted for. 

4 • 8 Difficulties at the Leeward Plane 

In solving three-dimensional laminar boundary layer 
equations, difficulties were encountered in obtaining solutions 
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at the leeward symmetry plane -where a complex interaction 
between the inviscid flow and the boundary layer occurs, 

(Moore, ref. 24). Moore showed that the cross-flow momentum 
equation has a unique asymptotic solution at certain condition 
and non-uniqueness occurs otherwise. This indeterminacy was 
attributed by Moore to a lack of previous history of the 
fluid which enters the leeward region from around the cone. 
Solutions in this region were studied by a number of authors, 
Moore (ref. 25) ; Cheng (ref. 26); Vvedenskaya (ref. 27); 

Libby (refs. 28-29), and Dwyer (ref. 30). Murdock (ref. 31) 
obtained solutions in the leeward region which were not 
known previously and determined conditions for which the 
boundary layer in this region is independent of the out-of- 
plane flow. He investigated cases for which complete boundary 
layer solution does not exist and concluded that the boundary 
layer model has a defect in this case which results in dis- 
continuous derivatives in the leeward plane. This is the 
region where the boundary layer model breaks down and the 
continuity of mass is not satisfied. 

It was found in this study that the boundary layer 
solution at the leeward plane does not exist, and that the 
region within which convergence cannot be obtained grows with 
angle of attack. It is believed that the boundary layer 
equations are inadequate in this region which was demonstra- 
ted by the inability to produce solutions except at small 
and moderate angles of attack. This seems to be confirmed 
by the fact that the boundary layer thickness does not grow 
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as fast as was observed experimentally in the leeward region. 


4 . 9 Solution to the Problem 

The governing equations to be solved are (2.23) through 
(2.25). Implicit within this set of three equations, the 
quantities A^, A 2 , ... , A^g (see equations (Al) through 
(A19)) need to be specified. These nineteen quantities de- 
termine all the inviscid and geometric restrictions on the 

* 


flow. However, the quantity £ has yet to be determined. 

In order for the quantities Ag , Ag , and A^,. to remain finite 

* 

at the tip of the sharp cone, 5 must be chosen to be /£ 

(this could have been deduced from the usual form of the 
similarity transformation as is used for most compressible 
cone solutions) . Note that now the three governing equations 
(2.23) through (2.25) are similar in 5 for the case of con- 
stant entropy since there is no dependence of the inviscid _ 
quantities on £ and since there is no explicit dependence of 
the set A^, , ... , A^< 

along with terms such as 


A^ , A -j , A-^, and Aj g are zero 
30 


— , , and . For this set of 

cases, injection will not be considered so that H w will be 
zero. Therefore, the boundary condition on <f> as in equation 
(2.34) reduces to the following: 


3 <p 

n 


& - - ♦ 


w 


w 


(4.28) 
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The problem of obtaining the starting solution for the 
step-by-step integration in the redirection around the cone 
still remains. Note that a similarity solution exists at 
the point at which the transverse pressure gradient, and 
likewise the transverse velocity given by equation (4.27) 
is zero. For a circular cone this occurs at the windward 
and leeward streamlines only. The integration will advance 
in the positive redirection from the windward streamline 
(at n = 0) in the case of a circular cone. 

Thus, the entire set of equations has been specified, 
and the finite-difference scheme has been developed so that 
the solution can be integrated step-by-step around the body 
(Fig. 3.2) . 

4 . 10 Results 

Based on the analytical model presented, a computer 
program for solving the system of parabolic boundary-layer 
equations for a circular cone was developed. This program 
is a modification of the one developed by McGowan and Davis 
(ref. 5) . This program begins computations for constant 
entropy at first, computes the boundary-layer thickness at 
the second step, and switches to the program with entropy 
swallowing at the third step and approximates the boundary 
layer thickness at this step (which is not known a priori) 
by the boundary layer thickness at the preceding step. 

To be more accurate, after the completion of computations 
over the prescribed number of steps in the first sweep, a 
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se co nd sweep can bo initiated by returning to the beginning 
of the program and using the true boundary layer thickness 
at each station, respectively from the first sweep. The 
option for performing the second sweep has been provided. 

As will be shown later, the error resulting from the approx- 
imation of the boundary layer thickness by the value from 
the previous step is negligible. 

The computed results include: skin friction, heat 

transfer, boundary layer thickness, displacement thickness, 
Stanton number and the normal direction profiles of velocities, 
temperature, density, and shear parameter. The results of 
the computations will now be discussed in detail for particular 
cases presented in Table I. Comparison of results with ex- 
perimental data and results from other sources will be made. 

The computations were performed at conditions which cor- 
respond to the tests performed by Tracy (ref. 23) and are: 


M = 8 

to 

6 = 10°, = 25° 

y c Y C 

T = 84 . 2°R 
00 

T = .46 T 
w o 


Re = 4.04xl0 4 where L = 1 ft. 
Y = 1.4 
Pr « .738 

R = 1718 ft 2 /sec 2o R 
C = 198. 6°R 


The computed cases are tabulated in Table I. The computations 
were performed on an IBM 360-65 computer and required a . 
storage capacity of 400K. For the computed cases, 76 points 
in normal direction from c, = 1 to c = 0 and 61 points in 
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truns verse direction n were used. It was found that this 

combination of steps represented an optimum from the point of 

view of convergence, accuracy, and computing time. It took 

approximately 2.3 seconds to compute the solution at a grid 

point. The compile time was 4 min. The skin friction and 

heat transfer were normalized by values C,,. and q at zero 

angle of attack, respectively (Table II) . Longitudinal and 

transverse skin friction and heat transfer distributions in 

the transverse direction for <j> = 10°, and a = 8°, 10°, and 

c 

12° and various values of the Reynolds number are shown in 
Figures 4.6 to 4.15. The same parameters are shown for <i> 

= 25° and a = 12.5° in Figures 4.16 and 4.17. These figures 
show values for the constant entropy cases (no swallowing) 
and for entropy swallowing where inviscid data were interpo- 
lated at a distance from the wall equal to the thickness of 
the boundary layer, <5 . Heat transfer results for <t> = 10° L 

are compared with the tests by Tracy (ref. 23) * 

The longitudinal and transverse skin friction distribu- 
tions for 1> =10°, a = 8°, and Re = 4.04 x 10^ are shown in 

c 

Figure 4.6. The Reynolds number is moderate and the boundary 
layer has a significant thickness. As a result of entropy 
swallowing the values of the skin friction are higher than 
those for constant entropy over almost the entire periphery: 
about 2% at the windward and 4% at the leeward region. A 
similar trend is displayed by the heat transfer in Figure 4.7. 
Heat transfer shows good agreement with the test values by 
Tracy (ref. 23) except toward the leeward region and windward 
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region where the computed value is about 3% higher than the 

test results. Figures 4.8 and 4.9 show the same parameters 

6 

for higher Reynolds number. Re = 4.04 x 10 . Because of the 

thinner boundary layer, smaller entropy swallowing effects 

are noticeable as is expected. Figures 4.10 and 4.11 show 

o 4 

the computational results for a = 10 and Re = 4.04 x 10 . 

As before, there is a marked effect of entropy swallowing 
on and q^. A fair agreement of computed heat trans- 

fers is still apparent, although the differences are larger 
near the planes of symmetry and are about 4% at the windward 
side. The effect of entropy swallowing is reduced by a high- 
er Reynolds number, as is noticeable in Figures 4.12 and 4.13. 
Considerable difficulties -were encountered in obtaining solu- 
tions for a/ <j> c > 1. This is the limit of validity of the 
inviscid data which becomes progressively more inaccurate as 
a increases. Figure 4.4 indicates the scatter of the computed 
inviscid data for a = 12°, <j> = 10° as obtained by Jones 

(ref. 14) and the distribution of the flow Mach number in the 
cross-flow plane, M c> It is obvious that the inviscid values 
are inaccurate and the fact that M reaches a value of unity 
at 0 = 95° indicates a possibility of internal imbedded shocks 
(ref. 27). Results for a = 12° are shown in Figures 4.14 
and 4.15. Figures 4.16 and 4.17 show results for <J> c = 25° 
and a = 12.5°. Although the trends are still the same, the 
differences between results for entropy swallowing and no 
swallowing seem to increase. Figure 4.18 represents a com- 
pilation of heat transfer distributions for a 10° cone at 
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a - 8° and includes experimental data by Tracy, results of 
computations by Boericke (ref. 4) based on experimental pre- 
ssure, and pressure by Moretti for M a = 7.95, and the results 
by the present method for M o = 8. The results of this analy- 
sis for entropy swallowing at M a = 8 are higher than the 
results based on pressure by Moretti at = 7.95 as should 
be expected. A similar comparison of heat transfer results 
for a = 12° is shown in Figure 4.19. 

The results obtained in this analysis can be compared 
with the recent solutions by Mayne (ref. 12). The results 
for (j) = 10° and a = 8° (Fig. 4.6, 4.7) show good agreement 

of skin friction values and heat transfer except in the 
windward area where Mayne ’’s results show that longitudinal 
skin friction is identical with the classical value and heat 
transfer which is lower than the classical value. In the 
case for <}> c = 25° and a = 12.5° (Fig. 4.16, 4.17), there is 
qualitative agreement in the skin friction, except that the 
results by Mayne for transverse skin friction is becoming 
lower than the classical toward the leeward side, while the 
heat transfer results by Mayne are slightly higher near the 
windward side and become lower for 0 > 60° than the classical 
values. The heat transfer results by Mayne for these two 
cases seem to lack consistency as compared with the results 
of this analysis which are consistently higher or identical 
for decreasing a values and high Reynolds numbers with stream 
line swallowing for the range of computed cases. It is be- 
lieved that these differences are due to the fact that the 
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streamline swallowing effects in the windward plane have been 

neglected by Mayne . Mayne argues that the vorticity effects 

are negligible in this region and assumes the similarity 

solution at £ = 0 , q = 0 to be valid for £ > 0 and n = 0 . 

Although the variation of flow properties in the windward 

plane is small, they are finite. The effects in the adjacent 

meridian planes are increasing and result in different cross- 

3 0 

wise derivatives, — , which lead to higher values of viscous 

3 n 

parameters in this area. 

The circumferential boundary layer thickness distribution 
for a = 8 °, 10°, and 12° is shown in Figure 4.20. The rate 
increase of the computed values of 6 toward the leeward plane 
is not as high as of the experimental data by Tracy. The 
streamwise boundary layer thickness distribution in the mer- 
idian planes 0 = 0°, 90°, and 150° is shown in Figure 4.21. 

It follows the law: 


O “ _ • 

/Si; 

The variation of the entropy parameter p/p Y in the circumfer- 
ential direction is. shown in Figure 4.22. The entropy decrea- 
se is larger for increasing angle of attack a and for decreas- 
ing Reynolds number (larger 5), which is consistent. It can 
be concluded that the entropy layer grows with the increase 
of a and 6 . 

It has been observed that the computed results fluctuate 
somewhat after the program switches to the variable entropy 
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scheme as shown in Figure 4.23. This is due to the impact 

of the different inviscid data set used in the variable 

entropy part of the program, different parameters used in 

the finite difference schemes and computed in the constant 

entropy phase of the program. This is also consistent with 

a slightly different inviscid data set and due to nonsimilar- 

a 0 

ity of the present scheme (~ ^ 0) . These fluctuations 

3 C 

practically vanish after the third step. 

It was pointed out that the entropy swallowing effects 
are accounted for by estimating the inviscid data at the 
distance from the wall equal to 6 from the previous step. 

This approximation results in some error which is shown for 
a = 10° in Figure 4.24. The results from the second sweep 
(second set of computations based on true 6 from the first 
set of computations) are higher by .1% at the windward side, 
.3% and .5% at 90° for and C^, respectively and go to 3% 
and 1% for a ■ and in the leeward area. For most practical 
purposes this error is small, and in view of the fact that 
the leeward results are not precise, the first sweep can be 
considered as satisfactory. 


4.11 Conclusions 

The results of the present method for solving the laminar 
boundary layer equations for a circular cone with entropy 
swallowing have demonstrated the validity of the approach of 
this analysis for the representation of viscous flow with the 
effects of entropy swallowing. In particular, the computed 
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heat transfer results < 3 vv / < 2 wo are in good agreement with the 
experimental results by Tracy and with results produced by 
other methods (ref. 4,12). The effects of entropy swallowing 
result in values higher than the values for constant entropy 
and this trend is consistent within the range of the (j> c , a 
and Re values used in these computations. Larger differences 
occur near the leeward region and this is believed to be due 
to the inadequacy of the boundary-layer equations as demon- 
strated in: 

1. The inability to predict a sharp rise of boundary layer 
thickness in this region 

2. The inability to produce a solution near the leeward 
plane which is deteriorating with increase in angle 
of attack a . 

The inviscid flow variables by Jones {ref. 14) seem to be 
satisfactory for a/<J> c up to unity and become inaccurate for 
larger relative a values over most of the body. They become 
inadequate near the leeward plane where strong interaction 
effects may be important. 

The error resulting from approximating the boundary 
layer thickness <5 by a value form the previous step is 
practically negligible and single sweep computations are 


sufficient. 
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5. SOLUTIONS FOR BLUNTED CONES AT ANGLE OF ATTACK 

Nearly all real configurations for hypersonic vehicles 
or lifting bodies require a bluntness of finite radius in 
order to reduce the heat transfer. Since a blunted body is 
of considerable importance, the second part of this work is 
devoted to the blunted-cone configuration, in particular to 
spherically blunted circular cones at angle of attack. The 
method of solution is similar to the one used for a sharp 
cone; consequently the same computer program can be used with 
proper modifications. 

5.1 The Geometry of the Problem 

The analysis of the viscous boundary layer flow around 
a blunted cone is performed as follows. The spherical part 
with symmetrical flow is solved first (Fig. 5.1, region I). 
After the computations have been carried out a certain distance 

, v» 

beyond the body axis, a transfer of coordinates to the spher- 
ical part about the body axis with unsymmetrical flow takes 
place (region II) . At the sphere-cone junction, the program 
switches to the conical coordinate system (region III) and 
computations are continued in the same way as on a sharp cone. 

From the expression for the square of a line element in 
an orthogonal curvilinear system 


•? 7 7 2 2 2 2 

dJT = dc + h^ dr, + h 3 dc 


(5.01) 
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one can write for a spherical system (Fig. 5.2) 


2 2 2 2 2 2 
d£ = dr + r d<j> + (r sin<j>; d8 z 


(5.02) 


and identify the scale factors for* the surface coordinates 
as 


h^ = r , h 2 = r sim}> , 
S = <J> and n = 0 . 


(5.03 a-d) 


When nondimensionalized by the radius of the sphere r , the 
following relations hold for the surface of the sphere 


(5.04 a-d) 


5.2 Flow About a Sphere-Cone 

With the scale factors h^ and h 2 identified, the invis- 
cid coefficients are determined for the symmetrical region 
I with special attention given to the case of the stagnation 
point (<j> = 0) and the region outside the stagnation point 

U > 0) . 

The program for the flow about the sphere at an angle of 
attack is broken up into two parts (Fig. 5.1): 


h, = 


= 1 


3h, 

' d 

W 


= cos<j> 


h_ = 


sin<t> and 


3 h, 

d 

frT 


= 0 
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a) symmetrical flow about the wind axis with zero crossflow 

3 0 

velocity v and azimuth symmetry (— = 0), this is denoted 
e a n 

by region I. 

b) unsymmetrical flow about the body axis with non-zero 

3 0 

crossflow velocity v and crossflow derivatives (-—=■ ^ 0) , 

e 9 n 

this is denoted by region II. 

Thus, these two regions must be treated separately. One 
set of inviscid coefficients A^ is specialized to region I 
and is used to obtain solutions up to an angle <J> = 3a. 

During this phase the results in the plane of symmetry from 
point ISI to ISF (Fig. 5.3) are stored and will be used as 
starting solutions for the unsymmetrical region II. Using 
the geometrical relations between the wind and body systems, 
the results at a finite number of points from ISI to ISF 
are interpolated evenly over the entire n region at IWEND 
points and are used to obtain solutions at the ISF+1 stream- 
wise station and up to the sphere-cone junction CJ. 

The position of a point and the sphere in terms of the 
body coordinate system (x,r,0) is related to the coordinate 
x in the symmetrical wind system for a given angle of attack 
a by the following expression (Fig. 5.4) 

X, = R + (x - R ) cosa + r sina cos 0 (5.05) 

k n o n o 

Thus, any point on the sphere in the body system along line 
1-1 can be reached from the wind system by determining the 
coordinate X, . 
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Since a sphere has no preferred orientation the flow 
is symmetric in the wind system. The flow in the body system 
is unsymmetrical . Any point on the sphere and the correspond- 
ing flow variables at this point can be obtained from the 
symmetrical flow about the sphere in the wind system by ro- 
tation about the wind axis. 


5 . 3 Method of Solution 

The method of solution of the laminar boundary-layer 
equations for a sphere-cone is essentially the same as the 
one applied to a sharp cone. For solving the symmetrical 
sphere flow in region I, it was established (equations 5.03, 
5.04) that 

3h 2 

c = <p h i = 1 Ytr = cos * 

3h 2 

n = 0 h„ = sin<j> - — = 0 

Z dr) 

and for similarity to exist at the stagnation point (Appendix 
C) 

l = 1 (5.06) 

With these relations the inviscid coefficients A^ are special- 
ized to region I by realizing that the crossflow velocity 

V =0 and crossflow derivatives are = 0 and that the 
e .3 n 

streamwise derivatives exist, ^ 0 . In order to be able 
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to transfer to region II where unsymmetrical flow exist, 
the computations are carried out to a value of £ = 3a. The 
computed parameters in the leeward plane are stored and in- 
terpolated over the half-range of n in the plane normal to 
body axis and passing through the last computed point in 
the leeward plane ISF and are used as starting values in the 
region II. This region is handled the same way as region I, 

except that the crossflow velocity and crossflow derivatives 

90 

are not zero (V ^ 0, ~ j* 0) . After the sohere-cone junc- 

6 d rj 

tion has been reached, computations are performed in the 

conical coordinate system as for a sharp cone, except that 
* 

K = 1 . 

5 • 4 In viscid Data - 

The inviscid flow variables for the sphere-cone analysis 
were obtained by the method of characteristics by Rakich 
(ref. 32) . The inviscid input data for region I consists of 
the pressure distribution p = p(<J>), for f = 0 to ^ = .7986 
radians (Fig. 5.5). With stagnation values specified one 
obtains the entropy parameter K = P q /p o Y and total enthalpy 

p 

o y 

H = — — hp- and velocity u = u (<J>) 
o p v-1 J ee Y 

o' 

u e 4 (h o ' <P e Y ' 1 K)T ’(^x) (5.07) 

In the region II about the body axis, the pressure distrib- 
ution is generated from the symmetrical pressure p(4>) from 
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rcqi.on X in 7 meridian planes spaced evenly 30° apart by 
means of relations from section 5.2. These pressure curves 
are matched with specified pressure distributions in the 
same meridian planes respectively, beginning at the first 
location at the axial coordinate x IIo . The longitudinal 
velocity in the windward plane is generated using the re- 
lation 



(5.08) 


The longitudinal velocity in each cross-plane is obtained 
by the New ton- Rap son iteration method described in section 
4.6 and the crossflow velocity is obtained from the £-momen- 
tum equation 

(5.09) 


u 3u_ , 3u 

V e - tan<j> (— jj — + iTH ^ ^ 


P V 

r o o 


IE) 


1 3p e 

The relation — — is computed from equation (5 .8) . An 

v e ^ 

alternative method for computing the crossflow velocity is 

determined by the following method. Bernoulli's equation 

(4.09) is differentiated with respect to £ and the derviative 
3p - 

is introduced into the momentum equation (2.01) which is 


35 


specialized to inviscid flow. After some rearrangement one 
obtains the irrotationality condition on the surface 


(h 1 u 



a 


( Ve> 


0. 


3 

3 n 


(5.10) 
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Using a two-point backward finite difference expression 
3V e 

for — — — one arrives at 

o t. 


V 

e i+l 



(5.11) 


Equation (5.11) is much easier to apply than (5.09). Equa- 
tion (5.10) must be modified if the surface of the body is 
not a surface of constant entropy. 


5 . 5 Results 

Computations for spherically blunted cone were performed 
for the following conditions: 


M = 8 

oo 



T = 418 . 87°R 

°o 

T = .6 T 
w o 

R = l.OxlO 6 
e 

R =. 1 ft. 
o 


v*- 


The results for region I are shown in Figure 5.6 and include 
skin friction, heat transfer, pressure, velocity, u g , and 
longitudinal derivatives of pressure and velocity- Skin 
friction and heat transfer results for region II and III 
are shown in Figures 5.7 and 5.8. The matching of computed 
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values at transition from region I to II and from region II 
to III is good. The circumferential distribution of the 
skin friction and heat transfer parameters in two cross- 
planes at <p = 56° and <j> = 92° is shown in Figure 5.9. Long- 
itudinal pressure and velocity distributions in 7 meridian 
planes is shown in Figure 5.10, and longitudinal derivatives 
of pressure and velocity are shown in Figures 5.11 and 5.12. 

A typical circumferential distribution at <f> = 56° of p , 

9v e ^ 

u , v and ^ is shown in Figure 5.13. The above results 

0 0 o T1 

for Crocco variables are compared with viscous flow solutions 
for laminar flow on a blunted cone at the same conditions 
using Levy-Lees variables (ref . 33) in Figures 5.7 and 5.8. 
The same pressure was used in both cases, however, the in- 
viscid data were handled differently. This probably accounts 
for some differences in results. Otherwise, the results 
show fair agreement as far as magnitude and distribution of 
parameters are concerned. 

5.6 Conclusions 

The solutions for viscous laminar flow over a spherically 
blunted cone were obtained by solving the conservation equa- 
tions in Crocco variables. The results show fair agreement 
with similar computations based on Levy-Lees variables and 
match fairly well at transition points from region I to II 
(unsymmetrical flow in body system) and from region II to 
III (conical flow) . 

Considerable difficulties were encountered in handling 
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the inviscid data, pressure, velocity and their derivatives. 

In particular it was difficult to obtain smooth distributions 
of these curves. Lack of continuity of these curves causes 
scatter of results or renders a solution impossible. There- 
fore, a more accurate set of inviscid data would be desirable, 
and work needs to be done on determining the best way to 
handle the inviscid data. 
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TABLE I. CONDITIONS OF COMPUTED SHARP CONE CASES 


Case 

M 

CO 

♦c (o) 

T /T 
w' o 

R 

e 

a (o) 

1 

8 

10 

.46 

4.04xl0 4 

8 

2 

8 

10 

.46 

4 . 04xl0 6 

8 

3 

8 

10 

.46 

4.04xl0 4 

10 

4 

8 

10 

.46 

4.04x10 s 

10 

5 

8 

10 

.46 

4 . 04xl0 4 

12 

6 

8 

25 

.46 

4.04xl0 4 

12.5 

' 
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TAIil.K II. COMPUTED SHARP CONE PARAMETERS AT a - 0° 


Case 

M 

CO 

j 

M o) 

ot.(o) 

T /T 
vr o 

; 

Cr: /xRe 
C 

^w /(pu3) « /xRe 

1 

8 

10 

0 

.46 

0.39925 

-.05537 
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8 

25 

0 

.46 

.71229 

-.11582 






















Figure 1.2 Streamline swallowing on a snhero-cone 
at zero angle of attack 
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Figure 1.3 Streamline swallowing on a sharp cone 

at angle of attack 
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Figure 4.8 Skin friction distribution on a sharp circular cone 
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Figure 4. 9 Ileat transfer distribution on a sham oimn 
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Figure 4.10 S3<.in friction distribution on a sharp circular cone 
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Figure 4.12 Skin friction distribution on a sharp circular cone 
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Fiqure 4.13 Heat Transfer Distribution on a Sharp Circular Cone 
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Figure 4.15 . Skin friction distribution on a sharp circular cone 
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Figure 4.16 Skin friction distribution on a sharp circular cone 
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Figure 4.17 Heat transfer distribution on a sharp circular cone 





Figure 4.18 Heat transfer distribution on 10 sharp circular cone at 
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Figure 4.19 Heat transfer sdistribution on 10 sharp circular cone at 




Figure 4.20 Circumferential boundary layer thickness distribution on a sharp circular cone 
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Figure 4.22 Circumferential entropy distribution on a sharp circular cone for the case of 
entropy swallowing a ’* : distance o 





Figure 4.23 Convergence of computed parameters after transfer to entropy swallowing 
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Figure 4.24 Circumferential error distribution between results from 2nd 
1st sweep for a sharp circular cone 





Figure 5.1 Sphere-cone geometry 



Figure 5.2 Spherical coordinate system 








and II for sphere 




Figure 5.6' Skin friction, heat transfer and flov; variables about a sphere 
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Figure 5.7 Longitudinal skin friction distribution in meridian planes for a sphere cone 
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Figure 5.8 Longitudinal heat transfer distribution in meridian planes for 
sphere cone 








Fiqure 5.10 Longitudinal pressure and velocity distribution in meridian planes 
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Figure 5.13 Circumferential distribution of p , u, v and - — , — on 
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